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O*. Abstract 

Q-T 

We formulate the Polyakov-Nambu-Jona-Lasinio (PNJL) model in 0+1 dimensions. The thermo- 



dynamics captured by the partition function yields a bulk pressure, as well as quark susceptibilities 
versus temperature that are similar to the ones in 3+1 dimensions. Around the transition tern- 



> 

r^. ; perature the behavior in the pressure and quark susceptibilities follows from the interplay between 

OO ' 

the lowest Matsubara frequency and the Polyakov line. The reduction to the lowest Matsubara 



frequency yields a matrix Model. In the presence of the Polyakov line the UV part of the Dirac 



o 

OO 

spectrum features oscillations when close to the transition temperature. 

X 



I. INTRODUCTION 

There has been a large success in modeling the finite temperature behavior of QCD 
using the Nambu-Jona-Lasinio (NJL) model |l|, |2j. The NJL model is based on an effective 
Lagrangian of relativistic quarks interacting through a local and chirally symmetric four- 
point interaction. It was suggested that this model may serve as a good approximation to 
the low-lying chiral excitations of the QCD vacuum as well as the QCD thermodynamics 
below the transition temperature T c . Key in the NJL model is the spontaneous breaking of 
chiral symmetry and the emergence of a chiral constituent quark mass, which is generated 
through the interaction of quarks with the chiral condensate. 

The main drawback of the NJL model is that it does not include the properties of color 
confinement. This leads to the problem that the model contains the wrong degrees of 
freedom near the transition temperature T c . This has led to the development of extended 
NJL models which include some effects of confinement by introducing the Polyakov loop as 



a new classical field which couples to quar 



cs. These models are referred to as Polyakov-loop- 



a new classical neld wnicn couples to quarks, inese models are relerred to as rolyakov-loop- 
extended NJL (PNJL) models [3j, U, |5|, |g, |7|, [sj] . Many aspects of these models have been 
extensively investigated recently, including thermodynamics and phase structure for two |9| , 
and three [lOj flavor systems, finite isospin systems ll|, imaginary chemical potential [121 ]. 
mesonic modes [13J and studies related to the fermionic sign problem and incorporation of 
fluctuations 14j . 

In this work we recast the PNJL model into a simple effective Lagrangian in 0+1 di- 
mensions. Modifications to the thermodynamics and susceptibilities as compared to the 
four dimensional case are discussed. We show that the key features of the four dimensional 
physics across the transition temperature are captured by the interplay of one Matsubara 
frequency against the Polyakov line. The model with one Matsubara frequency reduces to a 
matrix model. The resulting Dirac spectrum oscillates near T c . This paper is organized as 
follows: in section 2 we formulate the model. In section 3 we derive the phase diagram, the 
bulk pressure and quark susceptibilities. In section 4 we detail the matrix model and de- 
rive the pertinent mean-field equations for the resolvent. The Dirac spectra are constructed 
numerically for temperatures across T c . Our conclusions are in section 5. 



II. MODEL 
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Motivated by the recent work on the PNJL model \4, U |& 13, l8|, |g, lift Oil, ll2, \M, \lM we 
consider the following schematic Lagrangian density in one-dimension, including a NJL four- 
fermion contact term and coupling to a constant temporal background gauge field, whose 
dynamics is encoded in the phenomenological potential U: 

U = ^(i^D 4 + tm + 2/i 74 )^ + 9 -{{^H) 2 + (^W) 2 ) + U(<f>[A], <f?[A],T) (1) 

In the above equation, ip a j are quark fields where a = 1,2, ...,iV c are color indices and 
/ = l,2,...,Nf are flavor indices. For simplicity Nf = 1 unless specified otherwise. D 4 = 
d T + iA 4 is the covariant deriviative, m is the bare quark mass and A 4 = GA 4 X a /2 where 
A 4 is the temporal component of the SU(3) gauge field, G is the gauge coupling, and A a are 
the Gell-Mann matrices. We consider scalar fermions in our work, therefore the 7 matrices 
in Eq. (1) are 2x2 matrices. The axial-anomaly and the effects of U(1)a breaking will be 
discussed elsewhere. 

The mean-field analysis of ([1]) is readily carried out by the bosonization procedure which 
consists in replacing the four-quark interaction with color-singlet auxiliary fields defined as 

P = -2ig 2 (^ L ), 
pt = -2ig 2 (fy R ), (2) 

so that 

£ 4 = ipRi(d 4 + /M + iA 4 )i) L + i) ] L i{d 4 + // + iA 4 )tp R 

+^(P f + m)i) L + i; R i(P + m)il> R + ^PP ] + U((/>[A], <j)*[A],T) , (3) 

in the chiral basis, ip = (i{)r,iPl)- The effective potential for the background gauge field is 
expressed in terms of the traced Polyakov loop. We work in the Polyakov gauge and take 
the gauge field A 4 as time- independent. The traced Polyakov loop is then expressed as 

In our gauge choice, the Polyakov loop matrix is diagonal and defined as 

i A 4 
L = exp[— ] = diag(e M , e iu \ e~ l{ui+U2) ) . (4) 



While L is gauge dependent, <p[A] is gauge invariant and in general complex valued. The 
potential for the Polyakov line U(<j),<j),T) satisfies the Z(3) center symmetry. At low tem- 
peratures we expect the potential to have a minimum at = 0. At temperatures above 
To it develops a minimum which gradually forces <p ~ y 1 as T — > oo. The potential in [7|, 
which is used here as well, was fit in order to reproduce the pure-gauge lattice data in 3+1 

dimensions: 

U(d>, 0, T) -b 2 (T) lA b 3 ^ 3 + - 3) + b, m (5) 



where 



T 2 rr 6 vr r ' 4 



6 2 (T) =a + a 1 ( |) +a 2 (^V + a 3 (§Y ■ ( 6 ) 



The coefficients were fit in |7j to the lattice data for pure gauge QCD thermodynamics. 
They are given as a® = 6.75, a± = —1.95, a 2 = 2.625, 03 = —7.44, 63 = 0.75, and 64 = 7.5. 
The partition function corresponding to the above action is given as 

Z = f[d^} [d^\ [d(j>] [#] [dP]e~ s , (7) 

where 



S = / C 4 d 4 x = V 3 £ 4 dr . (8) 

Making use of the anti-periodicity of the quark-fields ip(r + (3) = —ip{r) and the fact that at 
finite temperature the operator id\ is invertible with a discrete spectrum uo n = (2n + l)-7rT 
the path integration over the fermionic fields can be done resulting in the following form of 
the partition function 

z=/ m i#-'-'* f "«nndeJ ,{m+p) ^^- a '\ (9) 

J n=-oo a=l yUn + lfi- A 4 l(m + P ] ) J 

III. THERMODYNAMICS 

We now discuss the thermal properties of the above model in the mean field (saddle point) 

approximation. The thermodynamic potential associated with equation (Q is 

00 
Q = U{<f>, 0, T) + Y,N C P 2 - T J2 Tr c \n[f3 2 u 2 + (3 2 (uj n + i{fi + iA 4 )) 2 }, (10) 



where S = l/{2g 2 ) 1 The above series does not converge and a divergent, temperature 



independent piece needs to be subtracted (see for example [15|). Following this procedure, 



the sum over n can be done explicitly, which up to an overall constant yields 

Q = U((j), 4>, T) + T,N C P 2 -N c u-T Tr c ln[l + Le'^-^] - T Tr c ln[l + LV^ + ^] , (11) 

where u = \P + m\. Using the identity Trln(X) = lndet(X) we arrive at the following form 
for the effective potential: 

Q = U((j), <p, T) + Y>N C P 2 -N c u - T ln[l + 3^'^-^ + 30e" 2/3( ^-^ + e" 3 ^^] 

- Tln[l + 30e-^ +/x) + 30e" 2 ^ +/x) + e ~ w{u]+ll) }. (12) 

In the saddle-point approximation, the values of P, <ft and <f) that maximize the above poten- 
tial are found by numerically solving the coupled system of equations: 

dtt dn dtt 

The solution of these equations then provides mean field values which can be used in evalu- 
ating any thermodynamic quantities, such as the pressure p(T,fi) = —Q(T,fi). The goal is 
two-fold. First, one would like to see if the above simplified PNJL model in 0+1 dimensions 
can reproduce the bulk properties of the PNJL model in four dimensions. Secondly, one 
would like to further reduce this model to a matrix model in zero dimensions by including 
only a finite number of matsubara frequencies in the sum of equation (TTOl) . The chiral con- 
densate as a function of temperature and traced Polyakov loop is plotted in figure [Q The 
dotted curve shows the mean-field trajectory for <fi. 

In the left panel of figure [2] we show both the chiral condensate and Polyakov loop 
as functions of temperature at /i = 0. The solid lines are the results from the model 
in + 1 dimensions, obtained by minimizing the thermodynamic potential of eq. ([121) . 



Figure 2 in reference 16J shows the same quantities for the model in four dimensions. 
We find qualitatively the same behavior for both the condensate, Polyakov loop, as well 
as the two susceptibilities (d{ipijj)/dT and d<p/dT). Also shown in these figures are the 



results using only the two lowest matsubara modes in equation ( TTOl as dashed curves. For 



1 The value of S is chosen such that the constituent quark mass at zero temperature is P(T = /i = 0) = ^ 
300 MeV, a value consistent with NJL models in 3+1 dimensions. 
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FIG. 1: Scaled chiral condensate, (ipip)/(ipTp)T=o, as a function of T and 4>. The dotted red line 
shows a schematic trajectory when 4>{T) is taken from the mean field calculation. 

temperatures T > 100 MeV the sum over the first two frequencies is a good approximation 
to the infinite sum. The vanishing of the chiral condensate at low temperature for the 
truncated frequencies is due to the occurence of (3 = 1/T in the weight factor in §§§, which 
vanishes as (3 = 1/T — > oo. When all Matsubara modes are included, this vanishingly small 
weight factor is overcome by the determinant part with infinitly many modes, leading to a 
finite chiral condensate at zero temperature as it should. We will come back to this point 
in the random matrix reduction. 

We now show the results for the quark number susceptibilites using both eq. (|12p . where 
the explicit sum has been carried out over all matsubara frequencies, and eq. ([TO]) where the 
sum includes only the two lowest frequencies (±7rT). The coefficients are extracted up to 
eighth order by a fit to the scaled pressure 

^ = 5><r>(£)"- ("I 

n=0 



First we should compare our results with those of reference 17| where this exercise was 
carried out in 3+1 dimensions. The solid curves in figures [3J17I show the scaled pressure and 
first four susceptibilities for the PNJL model in 0+1 dimensions where the sum is performed 
over all frequencies. Qualitatively, similar behavior is obtained in both the 0+1 dimensional 
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FIG. 2: (Color online) Left: Scaled chiral condensate ('tpip) /(V ; V ; )r=o (red curve) and Polyakov 
loop (blue curve). Right: d(ipiJj}/dT (red) and dcft/dT (blue). In both figures, the solid curves are 
evaluated by using the sum over all matsubara frequencies and the dashed curves by using the two 
lowest frequencies, ±nT. 



model used in this work and the four dimensional model used in 17|. Near the transition 
temperature T c the peak structures are again similar in both models as seen in ci through 
Cs, which shows the direct interplay between the Polyakov line and the lowest Matsubara 
frequencies. There are qualitative differences in the high temperature behavior. This is due 
to differences in the dimensionality of the problem which we now discuss. 

In order to understand the effect of fewer dimensions and see if there are any qualitative 
differences between the model in four and one dimension we look at the high temperature 
limit (no longer mean field) of an ideal gas of quarks and anti-quarks. In four dimensions 
the pressure is given as 



P_ N c N f 

T 4 7T 2 T 3 



poo 

/ dkk 2 \n[(l-n)(l-n)] 
Jo 



(15) 



where n = 1/(1 + e -(v^?+ri/T) anc i n = 1/(1 + e -Wk*+m*-ri/T} which leadg to a finite 
series in chemical potential for massless quarks 



P 



7tt 2 



/i\4 



Q K T> 127i 2{ T } 



(16) 



N c N f T* 180 
Immediately one can extract the high temperature behavior for the susceptibilities in four 



dimensions: c<i 



N c N f 



C\ 



N c N f 
12vr 2 



and cq = eg = 0. However, in the 0+1 dimensional NJL 



model, the pressure is given as 

| = iV c iV71n[(l-n)(l-7i)], (17) 

which for massless quarks at high temperature (i.e. n,n = 1/[1 + e T/i / T ]) leads to the 
following result 

P = , 1//X\2 1 /^s4 1 //ix6 17 ,/lv 8 , . 

iV c JV/T n 4VT' 96 M 1440 VT' 322560 M l j 

This leads to different asymptotic susceptibilities in 0+1 dimensions: C2 — > c 4 / , and now 

C4 — > |g^ has changed sign. Also, c% and eg are non- vanishing. Therefore, at least in the 

high temperature limit, one should expect qualitative differences between the model in four 
and one dimensions. 

The dashed curves in figures [SK show the same result using only the two lowest Matsub- 
ara frequencies in the energy sum. At temperatures close to T c the qualitative structure of 
the susceptibilities is reproduced. At higher temperatures the finite sum result approaches 
the full result. We note that in the current model the high temperature limit of the suscep- 
tibilities is never reached. This is due to the fact that and are treated as independent 
variables and, as pointed out in [8] , this tends to overestimate the difference between and 
4>. If we set 4> = 4>, the ideal gas result would be obtained in the high temperature limit. 



IV. QUARK SPECTRUM 

The present model can be simplified further by putting the left and right handed quarks 
on a discrete grid spanned by the spatial variable x = 1, 2, ...N x and choosing the auxiliary 
field P to be a constant in space and time. In frequency space this sets the restriction that 
only certain Matsubara modes can interact (n = m,k = I). 



/ dr((^W + 0AW) 2 ) = 4/? V Sn^k^ljlRm^U. (19) 

Jo ■ , 



"P 

'0 . . 

n,m,k,l 

One can now bosonize quark pairs of opposite chirality in eq. ([3]) using the auxiliary matrix 
Wnm = ^RmYLn where the upper indices refer to three-space and the lower indices to 
frequency space resulting in the following Lagrangian 



H 
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FIG. 3: Scaled pressure normalized to the ideal gas result. The solid line is from eq. (|12p while the 
dashed line is the result keeping only the two lowest matsubara frequencies, ±irT. 



U = ^74+^+2/i74-A 4 74)^+£iVcTr x ^ 

(20) 
where Q = u n t n <8> Is ® H-c- The four- Fermi interaction causes the quarks to interact as if 
they were moving in a random Gaussian potential provided by the auxiliary fields 



P = -2ig 2 (^ L ) 



(21) 



Note that in the above Lagrangian we have defined E = /3E = Vi/(2g 2 ). This was done 
in order to make a connection with the standard RMM used in the literature. We will first 
look at this standard model in the thermodynamic limit where n = V^/N = const. Then 
we will relax this assumption and include the temperature dependence in the action (i.e. 
n = V3/N = const.). Note that the correct temperature weight is required in order for the 
matrix model to reproduce the mean field results with the lowest two matsubaras. 




FIG. 4: C2 as a function of T/T c normalized to the ideal gas value. 



0.8 



c? 0.4 




T 1 1 1 1 1 1 1 1 - 



a i 



2.0 



2.5 



FIG. 5: C4 as a function of T/T c 
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FIG. 6: cq as a function of T/T c 
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FIG. 7: cs as a function of T/T c . 



11 



We first set <f)[A] to a fixed value so the potential U will not affect the dynamics. The 
resulting form of the partition function allows for the investigation of the quark spectrum 
in the presence of the background gauge field. 

V[W]e- NSTrWWi T[det\ P . (22) 

jj[ \ iW ] + iu n + /i - L4 4 m J 

In the above model we have set the Dyson index to two corresponding to N c = 3 and the 
matrix elements correspond to the chiral unitary ensamble (%GUE). Each matrix W has 
N = 2N UJ N X N C entries whose interaction matrix elements Wij are drawn from a gaussian 
distribution having variance E = 1. 

Without the inclusion of the A± term in the above random matrix model (RMM), the 
above partition function is the chiral random matrix model of [l9|. The addition of the 
background gauge field A4 serves as an imaginary color chemical potential. A similar model 
was considered in 20j where a non-random component was added to the Dirac matrix in 
order to simulate the formation of instanton-anti-instanton pairs. 

We now examine the above matrix model for Nf = 1 and m = \x = and restrict the 
frequency space to the two lowest Matsubara modes. In this case (ft = 4> and the eigenvalues 
will be real. The matrix model is composed of a random part R and a deterministic part 
D. The model can be re- written as 

V[R]e~ Nf]Tl ^ RR 'd& x . r! . s Q m) 

where 

/ d\ / r\ 

Q = + (24) 

^Doy \ Rt °/ 

and D = \ x ® diag(7rT + u, nT, ttT — v, —ixT + u, —irT, —nT — z/)with v = T arccos(-^— ). 

In the mean-field approximation, the resolvent for the RMM follows readily from the 
use of Blue's functions (B) which are the inverse of Green's functions or resolvents (G), ie 
B{G) = G(B) = z 22:]. The Blue's function for the random part is Br = z + 1/z, while 
that of the deterministic part is 



B D {z) ( . 

71=1 



12 



stl^DT (25) 



where D n represents the n th diagonal entry of the matrix D. The Blue's function for the 
RMM follows from the self-energy addition rule -Br+d = Br + Br, — 1/z. The Green's 
function or resolvent for the RMM follows from the inverse rule Br + d(G) = z, 

which is a seventh order algebraic equation for G(z), 

G 7 + a 6 G 6 + a 5 G 5 + +a 4 G 4 + a 3 G 3 + a 2 G 2 + a x G + a = (27) 

with 

a 6 = —6z 

a 5 = 1 - 3tt 2 T 2 - 2z/ 2 + 15^ 2 

a 4 = z(-5 + 12vr 2 T 2 + 8z/ 2 - 20z 2 ) 

a 3 = 3vr 4 T 4 + z/ 4 + 5^(2 + 3z 2 ) - 2tt 2 T 2 (1 + 9^ 2 ) - 4/3z/ 2 (l + 9z 2 ) 

a 2 = -2^(3vr 4 T 4 + v A - 3vr 2 T 2 (l + 2z 2 ) - 2z/ 2 (l + 2^ 2 ) + z 2 (5 + 3z 2 )) 

0l = -7r 6 T 6 + z/ 4 (l/3 + z 2 )-2z/V(2 + ,2 2 ) + 2 4 (5 + ^) 

+vr 4 T 4 (l + 2z/ 2 + 3^ 2 ) - tt 2 T> 4 + 3z 2 (2 + z 2 )) 
a = -^/3(3tt 4 T 4 + z/ 4 - 6tt 2 TV - 4i/V + 3z 4 ) 

The spectral density follows from G(z) through its discontinuity along the real axis 



(28) 



p(A) = — lim ImG(A + ie) (29) 

71 e-s-0 

The algebraic solutions of ( J27l) leading to (|29l) will be discussed elsewhere. 

Now we discuss the case when the explicit temperature dependence is included in the 
Gaussian weight. The model is now written as 

Z= f V[R]e- NpsFt ™ NMt det x , n , N Q (30) 

where £ = V3/(2g 2 ). The prior result for the resolvent corresponding to the standard 
matrix model can be re-used for the above matrix model via the re-scaling of T — > yT and 

p{\) -> v^p(Av^)- 

13 



A key difference with the standard chiral RMM 21( is the fact that the Gaussian R- weight 
factor has an explicit fi = 1/T which causes the Gaussian weight to weaken as we cross the 
transition temperature from above going to low temperatures. It is this mechanicsm which 
caused the chiral condenstate to vanish at zero temperature in the mean field analysis of 
the two lowest matsubara modes. Our RMM is suited for studying spectra near T c as it 
embodies the extra suppression in T encoded in the time-integration. 

Before we show the spectra we must discuss what values of the background potential to 
use. The most physical choice of to take when computing spectra would be along the 
trajectory shown in figure [U Instead, however, we simply choose to show spectra using 
= 0, 1 in order to demonstrate the maximum effect of the Polyakov line. The spectra are 
shown in figure [S] as functions of temperature for = 1 (left figure) and = (right figure). 
The inclusion of the Polyakov line causes the spectra to split into separate domains as seen 
in the right figure. 




3 -20 



-20 



FIG. 8: Spectral function as a function of T for 0=1 (left) and = (right). 



Of course, all of the above values of T and are not realized in nature. For the case of zero 
chemical potential (which is what is examined in this work) the value of A 4 = diag(z/, 0, —v) 
where v = Tarccos(^p-). At low temperatures 0^0 and v ex T. At higher temperatures 
— > 1 and v — » as dictated by the arccos dependence. The maximum v occurs a little 
above Tc- We therefore expect to see the strongest changes in the eigenvalues near Tc- 

Finally, we compute the eigenvalue density by taking an ensemble average over 0, which is 
the analogue of integrating over the large gauge configurations across the transition temper- 
ature. This is done at three temperatures, 0.75T C ,T C and 1.5T C . The distribution function 
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for is found from the potential U ((/>). At T c the probability distribution function for <fi is 
peaked at <j) = and decreases monotonically to zero at <fi = 1. Below T c there is even more 
strength at lower values of <fi. Above T c the distribution develops a peak at a finite value of 
<f). The result of this ensemble averaging is shown in fig 




P(A) 




FIG. 9: Spectral density of Dirac eigenvalues generated from an ensemble of 2000 240 x 240 
matrices at T = 0.75Tc (top) , T = Tc (middle) and T = 1.5Tc (bottom). The value of <f> used in 
each matrix was sampled from the distribution given by U((p,T). 

We find that the inclusion of a background gauge field brings about qualitative differ- 
ences in the macroscopic spectral density in comparison to the standard chiral random 
matrix model. There are oscillations present in the bulk of the spectrum that increase with 
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temperature. These oscillations are most likely above the Thouless energy and therefore in 
the diffusive regime of the lattice data. However, there are also qualitative changes in the 
low energy eigenvalues. In the standard chiral random matrix model the spectral density is 
a semi-circle even above T c . Therefore the slope of the spectrum at low energy is infinite. 
The inclusion of the background field changes this picture. At T c one can see that the slope 
of the eigenvalue distribution at A = is finite. 

The Polyakov loop raises the critical temperature from the standard NJL model. The 
mechanism for this is now clear by looking at the spectrum. Near T c one makes an ensemble 
average over all values of <fi. When the gauge field is present, the spectrum is modified and 
the eigenvalues are shuffled to lower values of A. Due to the weighting from U((f>,T), only a 
small portion of the eigenvalues is shifted to lower A. This gives rise to the slope seen in the 
fully integrated spectrum. These additional low-lying eigenvalues therefore shift the critical 
temperature to higher values due to the Banks-Casher relation. 

V. CONCLUSIONS 

We have constructed a simplified version of the PNJL model in 0+1 dimensions that em- 
bodies the essentials of the model in 3+1 dimensions. The bulk pressure and susceptibilities 
across the transition temperature are shown to follow from the interplay between the lowest 
Matsubara mode and the Polyakov line. Many features of the quark number susceptibilities 
in the flavor symmetric case are analogous to the ones observed in full fledged lattice simu- 
lations, providing simple insights to the dynamics at works in QCD. The flavor asymmetric 
results both at finite temperature and chemical potential will be discussed in a forthcoming 
paper. 

The PNJL model in 0+1 dimensions yields a simple chiral RMM with a Polyakov line. At 
finite temperature, the presence of the Polyakov line causes the macroscopic Dirac spectrum 
to be multihumped across the transition temperature, which is a consequence of the Z(3) 
symmetry breaking at high temperature in QCD. The oscillations are smeared but not 
eliminated by averaging over the Z(3) vacuua, a feature that should be accessible to lattice 
simulations in the diffusive or chiral regime. The effects of Z(3) breaking on the spectra at 
finite chemical potential will be discussed elsewhere. 
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